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New technologies providing tight focusing lens and mirrors with large numerical apertures and electro-optic 
modulation of single photons are now available for the investigation of photon-atom interactions without a 
cavity. From the theoretical side, new models must be developed which take into account the intrinsic open 
system aspect of the problem as well as details about atomic relative positions and temporal (spectral) features 
of the pulses. In this paper, we investigate the generation of nonclassical correlations between two atoms in 
free space interacting with a quantized pulse containing just one single photon. For this purpose, we develop a 
general theoretical approach which allows us to find the equations of motion for single- and two-body atomic 
operators from which we obtained the system density operator. We show that one single photon is capable of 
promoting substantial entanglement between two initially ground state atoms, and we analyze the dependence 
of these correlations on atomic relative positions (vacuum effects) and spectral features of the single photon 
pulse. 

PACS numbers: 42.50.Ct, 42.50.Ex, 03.67.-a 

I. INTRODUCTION 

The strong coupling between photons and atoms has been demonstrated in many different experiments in the context of cavity 
quantum electrodynamics (CQED). In one typical physical setting, Rydberg atoms with large electric dipoles couple very effi- 
ciently to a single photon enclosed in a high quality cavity sustaining a resonant electromagnetic field mode [1]. These favorable 
conditions cause the atom to couple strongly to a particular mode of the quantized electromagnetic field in the cavity. In prin- 
ciple, this strong coupling between single photons and atoms could be used to implement distributed quantum communication 
(quantum networks) where localized information carriers (atoms) exchange information by means of flying qubits (photons) J2j]. 
However, the cavities used to coherently couple photons and atoms possess very long photon storage times [1] making photon 
extraction inefficient. Even in an intermediate coupling regime where photons may easily leak the cavity and information could 
in principle be transferred to another node (cavity with atoms) J2l, the scaling of the multi-cavity system is experimentally hard 
to achieve. 

This problem has been motivating the investigation of possible experimental settings to promote the strong coupling of a 
single photon with atoms in free space, avoiding the complications related to the use of cavities for quantum networking lU. 
On one hand, there is the question of what kind of field mode in free space would offer the highest electric-field intensity at the 
atom for an incoming beam. This question has already been answered |0] and it turns out to be an eletric-dipole wave, i.e., the 
time reversed version of the wave originally emitted by an oscillating electric dipole at the point of interest (an atom). From 
this fact, it follows that in order to achieve the strongest atom-electromagnetic field coupling in free space, one should build an 
optical system (mirrors, lens, etc.) capable of creating incoming light waves with great overlap with those dipole waves. The 
application of deep parabolic mirrors and other optical settings for such a task has been reported H. 

From a quantum mechanical point of view, where the concept of photon is meaningful, one should ask about the existence of 
an initial state of the free space single-photon pulse capable of perfectly driving an atom from its ground state to an excited state 
(a 7T— pulse). In fact, this state exists and its properties are intrinsically related to the dipole wave mentioned above in the scope 
of classical electromagnetism |6jj. The single photon state capable of perfectly exciting the atom in free space is the time reversed 
photon state originating from spontaneous decay of a two-level system, and its temporal shape is that of a rising exponential 1 6j] . 
From the experimental side, electro-optic modulation of single photons provide the key ingredient to study the transient response 
of atoms to different single-photon waveforms Q. All these recent developments in addition to the experimental demonstration 
of strong coupling between weak coherent light and a single trapped atom without a cavity iU make the investigation of models 
and effects related to single photon excitation of atoms in free space a very timely field of investigation. Following these lines, 
we present here a model to treat the interaction of two atoms with a quantized pulse in free space, and we study the power of a 
single photon to create entanglement between the atoms. It is important to remark that when two quantum subsystems interact 
with a commom bath, such as the vacuum electromagnetic field, they may become entangled depending strongly on their initial 
state [9]. In this sense, a classical laser field may be used to prepare a suitable initial state and then produce entanglement 
between the atoms. This is so because atom-atom interactions appear as a consequence of the free space vacuum environment, 
as seen previously in the context spontaneous emission of a collection of atoms in free space Il0l - [l2ll . This is a passive approach 
to induce entanglement in the system. One may also want to attempt an active approach where the dynamics of entanglement 
may be somehow externally controlled. For a two-atom system in the presence of vacuum as a commom reservoir, a classical 
laser field can be used to controllably populate collective states (symmetric and antisymmetric) and build up a considerable 
amount of entanglement between the atoms lfl3ll . In the present work, we extend this active control paradigm to an all quantized 
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approach which, for the single photon level, must be more accurate than the approach considering a semiclassical model with an 
attenuated laser beam. Besides, by building a fully quantized model, one can explore quantum resources not present in classical 
electromagnetic pulses, such as light squeezing or entanglement to interact with the atomic system in free space. 



II. MODEL 



The Hamiltonian for the system consisting of two two-level atoms sitting at positions ri and r2 and interacting with the 
quantized pulse can be written (in the interaction picture and RWA approximation) as 



H 1 = -ihJ2 J rfk[a kiA ( 5l;k , A a+ + . 92;k , A( r+)e- 1 ^-")* - H.c. 



(1) 



where erf = \e)i (g | with i = 1 , 2 are raising atomic operators. We consider here the case where the atomic transition frequencies 
of both atoms are equal and given by w. In the absence of any mode-selecting cavity, all modes of the electromagnetic radiation 
field couple to the atom, and the coupling constant c/.; ;kA is given by 



where e k A are unit polarization vectors, Uu,x( r ) are the spatial mode functions and fi { are electric dipole moments of the atoms. 

In order to evaluate the quantum correlations driven by a single photon we need to obtain the dynamical evolution of the 
atomic density operator p(t). We achieved that by we employing Heisenberg equations for treating this problem, generali zing 
the single atom approach discussed in Hill . In some respect, this approach can be seem as a form of input-output theory lll5ll 
which is largely known for providing tools for the treatment of general baths. The atomic density operator can be obtained from 
solutions of reduced Heisenberg equations for one- and two-body atomic operators having the field variables traced out. In the 
Schrodinger picture, the density operator for two two-level atoms can be written as lfl6ll 

1 3 
p(t) = -[I (8)7 + d(t) -o-i(8)/ + /® s(i)-<r 2 + O n , m (t)<r? <8 of 1 ], (3) 

n,m— 1 

where / stands for the identity operator for each two level atom, er^ = {erf , erf, erf} = {a} , erf , of}, and the coefficients are 
calculated as expectation values dj = Tr(pa\), s,; = Tr(pa|) and O nm = Tr(pcr™cr™). We may evaluate the expectation values 
using the Heisenberg picture and then plug the results in (|3) to obtain the density operator in Schrodinger picture. 

We then aim now to solve the relevant set of Heisenberg equations for the whole system (atoms +pulse) prepared in the initial 
state | "J) = .92, where the single photon state reads IU7I1 

IU = E / dk so*k,AiW4,Alo>, (4) 

A J 

and the atoms are initially in the ground state. The factor .go ; k.A is given by (0 and evaluated at the position r representing the 
focal point. This corresponds, for example, to the field state after spontaneous emission by an atom at ro or an incoming wave 
focused on an fictional atom at r . It all depends on the choice of f(oJk) which is linked to the temporal features of the pulse by 



dujkfMe- 1 ^-^. (5) 

We will treat here two representative single-photon waveforms that can be readily generated using electro-optic modulation 0], 
namely the rising exponential pulse and the Gaussian pulse. They respectively correspond to the choices 



and 



j- / 1 \ j v/ftcxp(ftt/2) for i<0 

* W ~ 1 for i>0 ' W 



£(*)=( ^) exp(-tt 2 i 2 /4). (7) 
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From now on we denote c/; ; k.A as g; but the dependence on polarization and wave vector must be properly taken into account 
in calculations using these quantities. With that initial state, the set of equations of motion for reduced Heisenberg operators for 
our problem reads (see Appendix lAV 

-71 (K) + 1) - 7ia(of of ) + A 12 (oX) - 2G„ (*)£«,(*)> 



err = 



dV 

±( 
dt 

d 

777 ( 



dt 

= -72(K) + 1) - 7iaK<f> - A 12 (af <t|) - 2G 1 (t)fc 1 (*). 
W x 2) = + TiXof of ) + ^(K) + (crf» + 7l2 (cr?af) + G (t)Go(t) + G 1 (t)Q 1 (t), 

fo?) = -i( 7 i + 72)K^) - ^(K) - (af)) + iG (t)V (t) - iG^V^t), 

M) = -(7i + 72)K^> - 7iK> - 72K) + 2 7 i 2 (crfcrf) - 2[G (t)Qi(t) + Gi(t)Q a (t)], 



(8) 



where JC q (t), T> q (i), and Q ? (t) are auxiliary functions defined in (1A131 in the appendixlAl A close look at (lAlOt in the appendix 
lAl allows one to see that (erf erf ) = (ofcr^) and (erf erf) = — (erf erf). The constants appearing in © are given by 



2tt 



A 



(9) 



1' 



= 2lT ^2 J dkgig*5(u k - w), i 7^ j, 



(10) 
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The delta functions here express the fact that <jj is essentially constant for frequencies of interest which are centered around the 
atomic transition frequency. The decay rates 7 j are readily evaluated and the results are the well known free space atomic decay 

constants 7, = ^ (^) where /Lt, = The Lamb-shifts that would appear as a consequence of Cauchy principal values 
in 7, have been discarded since their contribution turned out to be insignificant compared to 7i , 7^ , and A 12 . However, this is 
not the case for the cross terms, where the principal value gives rise to a free space field-induced dipole-dipole coupling with 
magnitude A12 ( fTTT i which is comparable to 7* and 7^. This induced coupling has already been predicted for the spontaneous 
emission of originally non interacting atoms IU0l - [l2ll . To calculate integrals ( TTOb and ( fTTT i. we set our coordinate system such 
that the direction of the z axis coincides with the direction the atomic separation vector = r% — r 3 - ITltl . and we set the atomic 
dipoles parallel and polarized in the xz plane as = fay) [sin a, 0, cos a) , where a is the angle between the dipole moments 
and Vij . The result is given by 



7ij 



(1 -3 cos 2 a) 



cosjSjj) _ sin(<%) 
(^) 2 (<%) 3 



(12) 



and 



A12 = ^V7i72 sin 2 a 




1 j 



^77177(1 - 3 cos a) i ^ 



cos(cy,) 

S 3 - 
'j 



5 3 - 

>j 



•Fi(<5, 



(13) 



where S l3 = 2irr l3 /\, Ji(0) = sin(0)Ci(0) - cos(0)Si(0), and J" 2 (0) = - sin(0)Si(0) - cos(0)Ci(0), with Si(0) [Ci(0)] sine 
[cosine] integral functions and = |ry |. Here A is the wavelength associated with the atomic transition of angular frequency 



Finally, the temporal features of the pulse appear in (|8) and (IA 1 3t through the functions Go (t) and Gi (t) which have been 
defined as 



G (t) 



(14) 



/7o V^o 

where is related to the spectral function f(u>k) in the initial field state (|4]i through the transform (0. Now we have all the 
tools necessary to analyze the physical properties of the atomic system under single photon excitation in free space, including 
the single-photon induced build up of quantum correlations between the atoms in free space. 



4 



III. QUANTUM CORRELATIONS 



Quantum mechanics allows the existence of correlations providing quantum information resources not present in classical 
systems. One example is quantum entanglement which is defined as the type of correlation that can not be created using local 
operations and classical communications. We aim now to study the arising of this kind of quantum correlation between the atoms 
due to their interaction with a single photon in free space. Since we are dealing with an open system, and the atomic state is 
certainly mixed, we will consider here the entanglement of formation which is a suitable entanglement measure for such cases. 
For two two-level systems, this entanglement measure is given by Hill 

EoF = h[(l + Vl-C 2 )/2], (15) 

where h(x) = ~x log 2 x — (1 — x) loe^fl — x) and C is the concurrence. The solution of (0 led to an atomic density operator 
which falls into the class of X states IU9I1 for which concurrence and then entanglement of formation can be readily obtained. 
From the density operator matrix elements in computational basis, concurrence of an X state can be calculated as 

C = 2max{0, |p 23 (*)| - ^ P n(t)p Ai (t), \p 1A (t)\ - ^W)}. (16) 

In the following plots, we will keep a = ir/2, 71 = 72 = 7 and ro at the midpoint between r\ and r 2 in order not to favor 
any of the two atoms. The initial system state in t — > —00 was \ty) = |<?i, #2, lp), as explained before. Figure[T]shows a typical 
dynamics for the entanglement between the atoms in interaction with a single photon in a rising exponential pulse or a Gaussian 
pulse. Each choice of fi/7 corresponds to an optimum value regarding the peak of the entanglement curve for a fixed atomic 
separation. We can see that the rising exponential pulse is able to create more entanglement between the atoms than the Gaussian 
pulse. This might be related to the fact that rising exponentials correspond to the reversed dipole waves which perfectly excites 
a single atom, but the analogy is not complete here due to the presence of a second atom which couples to the first through the 
quantized electromagnetic field. Also, by comparing (a) and (b) in figure Q] one can notice the great sensitivity to variation of 
pulse bandwidths il in © and (0. The whole atomic response to the single photon pulse must also change strongly depending 
on the atomic relative separation r 12 through the collective damping rate (TTZl i and free space field-induced coupling constant 
dol l. Figure |2]illustrates this point by showing the maximally achieved entanglement (peak value) for several atomic separations 
ri 2 /A as a function of the pulse bandwidth O/7. One can see from figure [2]that it is indeed possible to optimize il in order to 
maximize the peak value of entanglement for a fixed relative distance ri 2 between the atoms, just like done in FigureQ] 




Figure 1 : Entanglement dynamics for initially unexcited atoms in interaction with a single photon pulse for different atomic separations and 
pulse shapes. In (a) the parameters used for the Gaussian pulse were (ru/X — 0.2, Q/'y = 2.55) and for the rising exponential pulse 
(?"12/A = 0.2, $1/7 = 1.77). In (b) the parameters used for the Gaussian and rising pulse were (ri2 = 0.2, O/7 = 1.0). 



Our calculation provides a quantitative and complete quantum electrodynamic description of the interaction between two two- 
level atoms and one single photon in free space. From this model it is possible to understand that the building up of entanglement 
between the atoms can not be solely explained using the qualitative picture which asserts that the atoms get entangled because 
they "see" the same photon in the absorption process. The situation is quite richer than that as exemplified in figure [2] Had 
the building up of entanglement been only dependent on this fact, the closer the atoms the bigger the entanglement would have 
become. However, plots in figure |2]clearly reveal that there is an optimum atomic separation for each pulse shape. For example, 
by fixing fi/7 w 4.0, one can see that r 12 /A = 0.2 is more favorable to entanglement generation (maximum entanglement) than 
f 12/ A = 0.1. This is a consequence of the fact that the atomic coupling ( fl~3T > and the collective spontaneous emission ( fT2l depend 
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Figure 2: Maximal entanglement for several atomic distances as a function of the normalized pulse bandwidth fi/7. (a) Rising exponential 
pulse (b) Gaussian pulse. 



on the atomic distance in an oscillatory manner, and the equations of movement depend on the pulse shape. Consequently there 
is a competition among several factors which determines the atomic response. All these factors can be completely understood 
only within a full quantum electrodynamic calculation such as the one presented here. 

It is also important to contrast our findings and previous work on the generation of atomic entanglement by detection of 
indistinguishable single-photon emissions I20I1 . especially the work by Cabrillo et al I2H1 . In this case, entanglement between 
the atoms is generated by interference effects and state projection accompanying a measurement. State projection is an essential 
ingredient in the scheme. In our study, there is no measurement whatsoever and our main goal is to show an alternative way 
to create entanglement among atoms which is due to the electromagnetic pulse and vacuum induced couplings. In our scheme 
the atoms are originally in the ground state while in the approach showed in ll2lll the atoms are initially prepared in the excited 
state by two classical laser fields. Consequently, there is a chance of measuring both atoms in excited state and this clearly 
distinguishing our proposal from theirs since here this joint probability is always zero. This follows from the fact that there is 
just one single photon in the driving quantum pulse. It may be argued that our proposal is not quite robust since spontaneous 
emission eventually destroys all the entanglement. However, it is precisely the emission of a photon by the entangled pair that 
can be used to propagate information between quantum nodes. The idea is to control the pulses optimally in the time window 
where correlations can be transferred from distant nodes by photons. This idea o quantum memory in free space is under current 
investigation ll22ll . 

It is worthwhile to compare the methods normally used in CQED to the ones available in these free space quantum optical 
setups. In CQED, atomic states are usually engineered by suitably choosing the initial state and the time of interaction with 
monochromatic fields. Here, the idea is to control the spectral properties of the pulse in interaction with atoms initially prepared 
in some fiducial initial state. We think that the results presented here may serve as a benchmark for new investigations about 
quantum state engineer and quantum communication in free space by using single photons as a tool. For example, it would be 
interesting to see whether or not a general state \<f>) = a\e\, #2) — /3 1 <?i , e-i) may be generated from the single photon pulse of 
equation (|4]i by suitably choosing f(u>k)- We believe that a generalization of the time reversal argument used in |6] to find the 
pulse state leading to a perfect 7r— pulse would allow one to solve this problem, but we have not tested this conjecture. 



IV. CONCLUSIONS 



We presented a fully quantized model for describing the interaction of single photons with atoms in free space. We showed 
that the free space correlated decay and dipole-dipole coupling between the atoms in conjunction with the single photon pulse 
can induce atomic quantum correlations. All features displayed by the atomic entanglement are strongly dependent on the atomic 
relative positions and the pulse profile, so that it is possible to achieve high control of the entanglement dynamics by carefully 
setting the system and pulse parameters. We believe our results can be important for quantum communication where the build 
up and distribution of atomic correlations by photon pulses is an essential ingredient. One possible future route is the use of 
a second pulse to redistribute the entanglement generated by the first pulse as predicted here. This might follow from results 
started in I22I1 . where the idea of quantum memory with a single atom in a half cavity is proposed. 

Our proposal has been made in the context of photonic qubits in free space but it can also be relevant to circuit quantum 
electrodynamics where strong extinction of the transmitted power in a ID transmission line coupled to an artificial atom has 
recently been achieved l23ll . In fact, cooperative effects, just like the ones presented here in the context of atomic physics, are 
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under theoretical and experimental investigation in circuit QED systems 12411 . and our results may find applications in these 
systems as well. 
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Appendix A: Equations of Motion 

In this appendix we present details about the derivation of the system of equations dS). From formal integration of the 
Heisenberg equation for a k x (t) obtained with Hamiltonian ([]]) one obtains 



(Al) 



On the other hand, a general atomic operator O evolving according to Hamiltonian ([TJ will obey 

J t O(t)=J2 J dk[(. 9l [a+(t),0(t)]+. 92 [a 2 + (t),0(t)]K iA (t) e -^"^ 



(A2) 



Even though atomic and field operators commute when taken at the same time, the choice of order here may simplify terms in 
further developments lEBTl . We found it convenient to write the operator ak,A(i) on the right and the operator a k x (t) on the 
left just like used in another context in 12511 on page 394. By substituting dAU in (IA21 > and changing the integration variable to 
t = t — t', one gets 



dt 



°(*)=EE / dk U[^(*):OW]«k,A(*o)e- i ^-^ i -a] £iX (ioK[^(0,0(t)] e +^--)* 



i=l A 

EE/ dk 

i,j=l A 



g* 9j [a+{t),0{t)] 



t-t 



dra+(t - T)e^-^g ig ^J{f), 0(t)} 



(A3) 



From now on, one must take care in order to not change the order in which atomic and field operators appear since now they 
appear evaluated at different times and then they do not commute. The last two terms in (IA31 > appear in the form of a integral 
over r of er ? + ^ \t — t) multiplied by the function 



/(*) = W dkgfae 



(A4) 



or its conjugate. The frequencies Wfe of the field modes cover a very wide range and gig* varies slowly with respect to w^. 
It follows that the set of oscillating exponentials of f(t) interfere destructively when r increases starting from 0, so that f(t) 

becomes negligible as soon as t > r c , where t c is the correlation time of the vacuum fluctuations. In addition, er ? (i — t) 
varies much more slowly with r, over a time scale Tr 3> r c , where Tr is the damping time of the atoms. We can then neglect 
the variation with r of (t — r) as compared that of f(t) in the integral ( IA3b . and replace ( 5 (t - t) by cr,+ ( } {t) which 
may then be removed from the integral. Assuming t — to ^ t c , we get in this way 



i=l A 



EE/ dk[g ig * X ^la+,0]a-(t) + g*g lX Ma+(t)[0,ar] 

i,j=l A J 



(A5) 
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with 



X (±) = 7T<5(w fe -w)±i 



. V 



LO k -U 



(A6) 



where V indicate the Cauchy Principal Value. In order to obtain dS), one needs to get (lA51 l projected onto the initial state 

IV'i) = l.9i,52, 1 P >, where 



|!p)=E / rfk 9o * k . A /K)< A (to)|0). 



(A7) 



All calculations are then performed in the Weisskopf-Wigner approximation where it is assumed that the coupling constant gi is 
constant for frequencies of interest centered around the atomic transition frequency u such that 



E J dkgia^x(t )e~ 



i{u k —iti)t 



|l p ) = Gi_i(t)|0), 



where 



7o 

and 7io, 7o an d €(t) were defined in (flQb . Q and (0. If O is Hermitian, we have from (IA5t that 

2 2 



(ft 



i = l 



=^2R G [G^ 1 (Vi|k + ,O]|^ )J + E E / dk2Re[^;x ( " ) <V'i|[^,0](77(t)|^ 1 ) 

i,i=l A ^ 



(A8) 



(A9) 



(A10) 



from which (0 can be readily obtained. The auxiliary functions fC q (t), T> q (t), and Q q (t) present in ^ have been defined as 

AC,(t)=2Re(^|o-^ 1 |Vi) (9 = 0,1), 
D (t) =2ilm(ip \afa2\ipi), 

Q (t) = 2Re(^okf^"l^i) ! 

Ci(*)=2Re(^okroal^i). (All) 
with IV^i) = |ffi,52jlp) an d IV'o) = |ffi,52,0 p ). The equations of motion for the auxiliary system can be obtained from 

,2 2 

-^o|0|Vi}=E[ G -i(^l^'°]W+ E E / rfk 5i5 * x ( -^V'o|[^,0]^(i)|V'o)+3^X (+) (V'o|<W[0 ! ^]|V'o) 

„' 1 A A 1 \ *J 

(A12) 



which follows from (IA5t and (lA8b . It results in the system 



/C g (f) 

e,(*) 



i 



.Ai5 



27^(*) + 2^2^(*)-^2 9 
1 



A 



-7 9 M*) + -7i2 Q 9 W - *-|^X> 9 - 2G g 

"(^7gffil +7?)^«(*) -7l2^gffil(*) ~ \p\l3 q (t) - lqJq®t(t) ~ -Tp^S 
"(^7gel +7g)Qg(*) -7l2Qgffil(*) ~ ^7l2^g(*) - 7g^gel(*) ~ —^^1 + 2G lQ 



(A13) 



with corresponding to addition modulo 2 (q assumes values in {0, 1}) and J q {t) = 2i luv{tp \a q ji/'i) 
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